Optical asymptotics via Weniger transformation 
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Starting from the resurgence equation discovered by Berry and Howls [M. V. Berry and C. Howls, 
"Hyperasymptotics for integrals with saddles," Proc. R. Soc. Lond. A 434, 657-675 (1991)], the 
Weniger transformation is here proposed as a natural, efficient, and straightforwardly implementable 
scheme for the efficient asymptotics evaluation of a class of integrals occurring in several areas 
of physics and, in particular, of optics. Preliminary numerical tests, carried out on the Pearcey 
function, provide a direct comparison between the performances of Weniger transformation and 
those of Hyperasymptotics, which seems to corroborate the theoretical predictions. We believe that 
Weniger transformation would be a very useful computational tool for the asymptotic treatment of 
several optical problems. 
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Physical asymptotics finds a natural field of applica- 
tions in optics. It is well known, in fact, that wave 
optics usually reduces to geometrical optics when the 
wavelength of the radiation is let to vanish, and that 
the result of such limiting procedure turns out to be sin- 
gular at special loci, called caustics, where unphysical 
intensity divergences are predicted^ The above transi- 
tion leads to a mathematical description of the wavefield 
given in terms of series expansions which turn out to 
be divergent^ Despite the opinion of the great mathe- 
matician Abel, 3 divergent series are fashinating objects 
which are considered tools of invaluable importance in 
physics. Within the asymptotics realm, diverging se- 
ries naturally arise when the highly oscillating integrals 
representing the wavefields are evaluated via standard 
asympotic methods, like the steepest-descent/i While it 
is a well known fact that the leading asymptotic expan- 
sion is given by the sum of all the contributions coming 
from the different saddle points belonging to the integra- 
tion pathji it is less known that all higher-orders asymp- 
totic corrections can be sistematicly evaluated and ar- 
ranged in the form of a diverging series^ Berry and Howls 
(BH henceforth}^ developed an interpretation scheme of 
the above divergence, termed hyperasymptotics (H hence- 
forth) , in which the divergent character of the asymptotic 
series obtained from the Dingle's procedure is regarded 
as a source of precious information, to be suitably de- 
coded, rather than a mathematical pathological behavior. 
The central result of BH is that such a decoding can be 
achieved by taking into account an intricate sequence of 
truncated asymptotic series involving the whole network 
of saddles, including also those which do not belong to 
the integration path. Since its introduction, H has been 
estensively used, up to recent timei^ 

The aim of the present Letter is to show that the 
decoding operated by H can also be achieved by using 
a nonlinear resummation scheme, the so-called Weniger 
transformation (WT henceforth)^ when is let acting di- 
rectly on the diverging series expansions built up only 



at the saddle points belonging at the steepest descent 
path. In particular, our idea of employing WT is based 
on the universal "factorial divided by power" law, ruling 
the asymptotic behavior of the Dingle's expanding coef- 
ficients, which naturally arises from the H scheme, and 
the fact that the same type of divergence is displayed by 
the Euler series, for which WT has been proved to be 
the most efficient resummation schemed We also present 
a direct numerical comparison between the results ob- 
tained by H and by WT when the wavefield associated 
with the so-called cusp diffraction catastrophe, 10 mathe- 
matically described by the Pearcey function, is asympot- 
ically evaluated. 

First of all, we briefly resume the main results of the 
Dingle's work and of its refinements obtained by BH. For 
simplicity we shall limit ourselves to the following class 
of phase integrals: 

1 = J exp[i$(s)]ds, (1) 

where <&{s) is a function of the variable s and C denotes a 
steepest-descent integration path, in the complex plane of 
s, passing through one or more saddle points, denoted by 
{s n }, which are solutions of the equation d$(s)/ds = 0. 
The integral in Eq. (QJ is then customarily written as the 
sum of separate contributions, say 2k, corresponding to 
the saddle, belonging to the integration path, formally 
expressed as^£ 

oc 

I*=exp(i$»)E T r W , ( 2 ) 

r=0 

where = $(sfc) and 

with Tfc denoting a small positive loop around the saddle 
Sfe. The fundamental result found by BH was that, al- 
though in principle only the saddle Sk does contribute to 
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the integral in Eq. ([2]), all other saddles Sh, with h 7^ k, 
thus including those not belonging to the path C, are in- 
volved through the following, formally exact, resurgence 
formula^ 

1 all saddles oo , , . 

h=l t=0 kh 

h^k 

where the (complex) quantities F^n = &k — termed 
singulants, were introduced, and 7^ 6 {0, 1} are bi- 
nary quantities, related to the topology of the saddles 
distribution^ The resurgence effect in Eq. Q constituted 
the basis for the development of H, and the same equation 
will be the starting point of our analysis. Equation (j4|) 

gives the asymptotic law of T$ for large r, i.e.,— 

T«ocM (5) 

where h* indicates the saddle corresponding to the min- 
imum value of |-Ffe/i*|. In particular, Eq. (O explains 
why the adiancent saddles are responsible of the diver- 
gent character of the series in Eq. {2}. The "factorial 
divided by power" asymptotic law in Eq. (|5|) is also re- 

(k) 

lated to the fact that the coefficients Tr are basically 
generated via a local expansion involving high deriva- 
tives at the saddle Sk, a feature which is common to all 
asymptotic methods.— ^ In particular, Eq. |5]) shows that 
the divergent character of the Dingle' series in Eq. ^ is 
equivalent to that of the so-called Euler series, defined 
as^ 

r=0 J ° 

which, on the other hand, is strictly related to the func- 
tional form of some quantities, called hyperterminants, 
which arc of pivotal importance in H£ Weniger showed 
that the Euler series in Eq. © can be summed directly 
from the sequence of its partial sums. 9 What is more im- 
portant, WT has proved to be the most efficient way for 
the resumming to be achived when compared to other 
resummation schemes, like for instance those involving 
Pade approximants. 9,13 We do not wish entering into the 
mathematical details about WT, for which we suggest 
the interested reader to consult the extensive literature 
(see for instance the recent review in Ref . [HI and Ref . [l5h . 
Our aim is rather to give a first numerical evidence about 
the retrieving capabilities of WT when is applied for the 
asymptotic evaluation of a class of integrals of fundamen- 
tal importance in optics, which have been treated with 
the use of H. This will also deserve to provide a numerical 
comparison between H and WT. 

Basically, WT operates a transformation on the se- 
quence of the partial sums of a diverging series to convert 
it into a new sequence which should converge to the gen- 
eralized limit of the series itself. In our case WT acts on 
the sequence, {S^}, of the partial sums of the series in 



Eq. J2]), and converts it into a new sequence, say {S m }, 
defined by 9 - 




The form of WT given in Eq. ([7]) can be roughly ex- 
plained as follows^ On supposing that our Dingle' series 
can be summed, the mth-order partials sum is written as 
the generalized limit (the "sum" of the series), say T«, 
plus a remainder, say Rm\ i.e., Sm = 7^ fc ) + R^\ for 
any index m. For a diverging series, is just the divergent 
character of the sequence {R$} that "hides" the finite 
limit T^ k f Thus the challenge consists in decoding the 
value of the limit from the exploding behavior of {Sm ^}- 
The strategy common to most nonlinear sequence trans- 
formations, like WT, consists in trying to compute an 
accurate approximation, say R^ , to the remainder RJ$ , 
by conceiving a suitable asymptotic model, in order for 

~(k) (k) 

R m to be eliminated from S m , which should provide a 
better approximation {sffl — Rm^} to the limit than 
{Sffl} itself. For the Euler series, and thus for several 
factorial divergent series characterized by the asymptotic 
law in Eq. ([5]) , such an optimal remainder model has been 
found^ and leads to the WT in Eq. Q. 

In the final part of the present Letter we present some 
numerical tests devoted to illustrate the performances 
of WT in comparison to H and also to put into evi- 
dence some limitations of WT when it is used in asymp- 
totics. In particular, we have chosen a class of phase 
integrals of pivotal importance in optics, namely the so- 
called cuspoid diffraction catastrophes , 10 for which the 
phase function is a nth-order polynomial of the form 
$(s) = s™ + YJk=i x kS k , with {xi,x 2 ,--,x n - 2 } being 
a vector of (possibly complex) parameters. Our example 
concerns with the asymptotics evaluation of the Pearcey 
function, corresponding to n — 4. Furthermore, we also 
let x\ = x and x 2 = y- The hypeasymptotic evaluation 
of the Pearcey function was used, as a preliminary nu- 
merical test, in the original paper by BH. 6 Quite recently, 
Davis and Kaminsk: 8 have also proposed an hyperasymp- 
totics computational scheme, based on Hadamard expan- 
sions, to retrieve very highly accurate evaluation of the 
Pearcey function for real values of x and y. There is 
also another reason for using such example to test the 
asymptotics of WT. Recently^ the two simplest exam- 
ples of diffraction catastrophes, the Airy (n = 3) and 
just the Pearcey functions, have been computed via WT, 
but starting from their ascending power series. Thus the 
results we are going to present will be useful for compar- 
ison between the use of WT as far as Taylor-like expan- 
sions and asymptotics expansions are concerned. The 
evaluation of the coefficients in Eq. §5§ can be exactly 
in closed-form, and it is reported, for the Pearcey func- 
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tion, in Ref. |6j in terms of Gegenbauer polynomials 
Figure Q] shows the behavior of the relative error achieved 

Relative Error 




WT Order 



FIG. 1: Behavior of the relative error achieved by the WT, 
versus the order m [see Eq. 0], when the Pearcey function is 
evaluated through the WT at the point (x — 7,y — \/2(l + i)) 
(dots) and at the point (x = — 7, y — 1) (solid curve). 

by the WT, versus the order m [see Eq. (0], when 
the Pearcey function is evaluated through WT at the 
point (x = 7, y = \/2(l + i)) (dots) and at the point 
(x = — 7,y = 1) (solid curve). The first case corresponds 
to that presented by BH in their original paper^ while 
the second one has been treated by Paris and Kaminski. 8 
The two cases differ by the number of saddles which con- 
tribute to the Pearcey integral, respectively one (com- 
plex) and three (reals). It is seen that in both cases the 
obtained accuracies are comparable to (and even better 
of) those obtained via the hyperasymptotics treatment. 
In particular, in the first case WT provides a relative 
error of the order of 10 -15 for m = 20, which is about 
5 degree of magnitude better of that provided, on equal 
terms, by H,~ 8 The ultimate relative error obtained via 



H (see Tab. I in Ref. Q) is of the order of 1CP 12 , obtained 
with m = 25, whereas WT provides, on equal terms, 
an error of 10~ 16 i.e., four degree of magnitude better. 
Moreover, in Ref. the same case was treated by using 
WT directly on the ascending power series of the Pearcey 
function, and it turned out that to achieve a relative error 
of 10~ 10 about 90 terms had to be considered in the par- 
tial sums sequence, whereas the same degree of accuracy 
is here reached with less than 10 terms. The example 
treated by Paris and Kaminski corresponds to the phys- 
ical situation of a point located inside the cusp-shaped 
caustics, where the intensity distribution associated to 
the wavefield displays a typical fringe pattern, and WT 
produces results which are as accurate as those find in 
Ref. la (see Tab. 1 of that paper). 

Before concluding the Letter, we want to give some 
warnings about possible limitations concerning the use 
of WT in asymptotics. We refer, in particular, to the 
fact that WT is not able, in the form given in Eq. 0, 
to resum nonalternating factorial divergent series, like 
the Euler series in Eq. ((6]) when z is real and positive*^ 
The above pathological situation occurs, for instance, 
when the Airy function Ai(z) is required at the so-called 
Stokes lines, corresponding to complex value of z such 
that |arg(z)| = 2tt/3,— where, on the contrary, the re- 
trieving capabilities of H basically remain unchanged. 
We thus expect, in the general case, that the effective- 
ness of WT will depend on the position of the complex 
vector \x\,x-2, ...,x n _2} with respect to the Stokes lines 
and surfaces! 20 ' 21 In this perspective, a deeper investiga- 
tion about the retrieving capabilities of WT is required, 
but the preliminary results here presented seem to cor- 
roborate our initial feeling that Weniger transformation 
would be of great usefulness for the asymptotic treatment 
of several optical problems. 
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